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We study the effect of a large magnetic field on the chiral and diquark condensates in a regime of 
moderately dense quark matter. Our focus is on the inter-dependence of the two condensates through 
non-perturbative quark mass and strong coupling effects, which we address in a 2-flavor Nambu- 
Jona-Lasinio (NJL) model. For magnetic fields eB < 0.01 GeV 2 (corresponding to B < 10 18 G), 
our results agree qualitatively with the zero-field study of Huang et al., who found a mixed broken 
^~>) ' ■ phase region where the chiral and superconducting gap are both non-zero. For eB > 0.01 GeV 2 

' and moderate diquark-to-scalar coupling ratio Gd/Gs, we find that the chiral and superconducting 

, transitions become weaker but with little change in either transition density. For large Gd/Gs 

however, such a large magnetic field disrupts the mixed broken phase region and changes a smooth 
^\ ' crossover found in the zero-field case to a first-order transition at neutron star interior densities. 
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I. INTRODUCTION 

■ The existence of deconfined quark matter in the dense interior of a neutron star is an interesting question that has 
spurred research in several new directions in nuclear astrophysics. On the theoretical side, it has been realized [l|,l 
that cold and dense quark matter must be in a superfluid state with many possible intervening phases 0, 0, S f 
between a few times nuclear matter density to asymptotically high density, where quarks and gluon interact weakly. 

■ The observational impact of these phases on neutron star properties can be varied and dramatic d, [Hi [H, 0, Q3 • 
^ ' Therefore, it is of interest to situate theoretical ideas and advances in our understanding of dense quark matter in the 

C*~) , context of neutron stars, which serve as unique astrophysical laboratories for such efforts. 

. The phenomenon of color superconductivity in Quantum Chromodynamics (QCD) at large quark chemical potential 

7— I ' fi is interesting for several reasons. Weakly attractive gluonic forces drive a BCS instability [15[ which results in the 

, formation of Cooper pairs of quarks, breaking several global symmetries spontaneously and resulting in the appearance 

' of (pseudo-)Nambu-Goldstone modes Quarks and gluons can acquire an energy gap whose magnitude depends on 

Q\ | the baryon density, pairing channel, quark masses, flavor of participating quarks, neutrality constraints etc (see [l6| for 

■ a recent comprehensive review). New massless modes arise in this phase, such as the superfluid phonon (corresponding 
£sl \ to gentle fluctuations of the baryon number) and a rotated photon (a linear combination of the vacuum U (l) em photon 

•*h , and gluonic hyperchargc). The collective response of such modes to external perturbations determines the transport 
/\ ' properties of the diquark phase, with many phenomenological consequences for neutron stars. 

jh ; 

At asymptotically high density /j, 3> Aqcd and for number of flavors Nf=3, the preferred pairing pattern is a flavor 
and color-democratic one termed the color-flavor- locked (CFL) phase ||. This idealized phase, while it displays the 
essentially novel features of the color superconducting state, is unlikely to apply to the bulk of neutron star matter, 
since even ten times nuclear matter saturation density (po) only corresponds to a quark chemical potential fi ~ 500 
MeV^ Ao cp (assuming 3 massless flavors). At these densities, quark mass effects can be important, and as pointed 
out in [17| . must be treated non-perturbatively. It is reasonable to think that the strange quark current mass, being 
much larger than that of the up and down quarks, inhibits pairing of strange quarks with light quarks. For the 
purpose of this work, we will therefore adopt the scenario of quark matter in the two-flavor superconducting phase, 
which breaks the color SU(3) symmetry to color SU(2), leaving light quarks of one color ( say "3") and all colors of 
the strange quark unpaired. Although this phase appears to be disfavored in compact stars [H, , we adopt it here 
to highlight the competition between the chiral and diquark condensates in the most straightforward way, without 
the additional complications of compact star constraints. Also, our results will be qualitatively true for the 2SC+S 
phase P. [20|. which can be studied similarly by simply embedding the strange quark, which is inert with respect to 
pairing, in the enlarged 3-fiavor space. 

Our objective in this letter is a numerical study of the competition between chiral and diquark condensates at 
moderately large \i and large magnetic field using the NJL model, similar in some respects to previous works [T3, US 
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HHj HH , which treat the quark mass non-perturbatively. Instanton-based calculations and random-matrix methods 
have also been employed in studying the interplay of condensates [H, [25|. In essence, smearing of the Fermi 
surface by diquark pairing can affect the onset of chiral symmetry restoration, which happens at fi ~ M q , where M q 
is the constituent quark mass scale J2(|. Since M q appears also in the (Nambu-Gorkov) quark propagators in the gap 
equations, a coupled analysis of chiral and diquark condensates is required. This was done for the 2-flavor case with 
a common chemical potential in (T3], but for zero magnetic field. We use a self-consistent approach to calculating 
the condensates from the coupled gap equations, and find small quantitative (but not qualitative) differences from 
the results of Huang et al [l7| for zero magnetic field. This small difference is most likely attributed to a difference 
in numerical procedures in solving the gap equations. We also address the physics of chiral and diquark condensates 
affected by large in-medium magnetic fields that are generated by circulating currents in the core of a neutron or 
hybrid star. Magnetic fields as large as 10 15 G are expected on the surface of magnetars and their interior field 
may be as large as 10 18 G, pushing the limits of structural stability of the star [2?], There is no Meissner effect 
for the rotated photon, which has only a small gluonic component, therefore, magnetic flux is hardly screened [29| . 
implying that studies of magnetic effects in color superconductivity are highly relevant [3(J EH- Including the magnetic 
interaction of the quarks with the external field leads to qualitatively different features in the competition between 
the two condensates, and this is the main result of our work. 

In SectionHH we state the model NJL Lagrangian and recast the partition function in terms of interpolating bosonic 
variables. In Section [ill! we obtain the gap equations for the chiral and diquark order parameters by minimizing the 
thermodynamic potential (we work at zero temperature throughout since T sta r <C fi). In Section llVl we discuss our 
numerical results for the coupled evolution of the condensates as a function of a single ratio of couplings, chemical 
potential and magnetic field before concluding in Section [V] 



II. LAGRANGIAN AND THERMODYNAMICS 



We employ a Lagrangian density for two quark flavors (Nf=2) applicable to the scalar-isoscalar and the pseudoscalar 
mesons and scalar diquarks 



C-NJL = £kin + £qq + £qq 5 £fcin = q(ift + M7° - ™>)q , 

£qq = G s [(qq) 2 + {qij5rq) 2 ] , C qq =G D (qij 5 efe c q c )(q c i-f 5 e f e c q), (1) 

where q = (u,d), /i is the common quark chemical potential 1 , m = diag(m„, rrid) is the current quark mass matrix 
in the flavor basis, f = (t l ,t 2 ,t 3 ) are the Pauli matrices in flavor space, (ef)ij and (e c ) Q/33 are antisymmetric 
matrices for flavor and color, while q° = —q T C and q c = Cq T with charge-conjugation matrix C = — i-f°~f 2 . Wc 
take m u — ma = mo 0. As the diquark coupling Gd and scalar coupling Gs depend on the form of the underlying 
interaction and are not universal, we fix Gs empirically and vary Gd/Gs as a parameter, within phenomenological 
bounds (2j|. Introducing the bosonic fields a : (qq), n : (qij^fq), A : (q c i j -f^€fe c q), and A* : (qij5€fe c q c ), the 
bosonized Lagrangian becomes 

Cb = q(ift + m° - m)q - q(<r + • r)q - ^A*(q c i~f 5 e f e c q) - ^(qij 5 e f e c q c )A - ° ^ ~ ■ (2) 

Chiral symmetry breaking and color superconductivity in the 2SC phase is manifest by non-vanishing VEVs for 
a and A (we do not include the possibility of pion condensation [Hj]). The partition function in the presence of an 
external magnetic field B is given by 

Z = N J [dgpg]exp{ J* dr J d 3 x(C b + £ em )} where C em = ^eQ(q4q - f M°) ~ ^ , (3) 
which can be rewritten in terms of the bosonized version as 



Our assumption of a common chemical potential is for simplicity; in an actual neutron star containing some fraction of neutral 2SC or 
2SC+S quark matter in beta-equilibrium, additional chemical potentials for electric charge and color hypercharges must be introduced 
in the NJL model. Furthermore there can be more than one diquark condensate and in general M u ^ ^ M s I2I.I . 
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Z c Z\/iZ,-$ where Z C ~N exp 



Zi, 2 = / [dQ][dQ}exp { / dr d 



[dg 3 ][dg 3 ]exp < J dr j d x 



dr / d A x 



a 2 A 2 
4Gl + 4G/7 + 8tt 



^kinCQ, Q c ) + \eQ{Q4Q - Q C 4Q C ) + ^QA~Q C + \q c ^ + Q 



(4) 



The bosonized part Z c serves as a constant multiplicative factor. The subscripts "1,2" refer to quarks of color 1 
and 2 with Q = q\% and "3" to quarks of color 3. The kinetic operators in Ckin{<l, Q c ) now read {%</) + pj° — M) where 
M = niQ + a. Also, (3 — 1/T, and we have used the notation A~(/A + ) = — i75e/e c A(/A*). Since the condensate 
of u and d quarks carries a net charge, there is a Meissncr effect for ordinary magnetism, while photon-gluon mixing 
leads to a massless photon. In flavorOcolor space in units of e — \figej \/Sg 2 + e 2 the rotated charge matrix is given 
by Q = Q ® I — / <g> T 8 /2-\/3 (T 3 plays no role; the degeneracy of colors 1 and 2 ensures that there is no long-range 
gluon-3 field). In our case, this translates to Q charges 1*1,2 = 1/2, di, 2 = —1/2, 1*3 = l,da = 0. With inert s-quarks, 
we also have si,2 = —1/2, S3 = 0. The gapped 2SC phase is Q-neutral, while overall charge neutrality of the matter 
requires a neutralizing background of strange quarks and/or electrons. In this article, we take the strange quark mass 
very large so that they do not play any dynamical role. We will limit ourselves to a discussion of the competition 
between condensates in a large magnetic field and not impose the charge neutrality and beta-equilibrium condition, 
which is known to stress the pairing and lead to gluon condensation and a strong gluo-magnetic field 

Evaluation of the partition function and the thermodynamic potential Q = —ThiZ/V is facilitated by introducing 
8-component Nambu-Gorkov spinors for each color and flavor of quark, leading to 



lnZ h2 = iln{Dct(/3G- 1 )}; lnZ 3 



llnlDet^Go 1 )}: 



(5) 



G~ 



0,Q J 

A+ 



0-Q 1 



0,Q J 





where [G^] 1 = (@ ± /17 + eQ — M) . The determinant computation is simplified by re-expressing the Q-charges 

in terms of charge projectors in the color-flavor basis, following techniques applied for the CFL phase [34[. With 
the color-flavor structure unraveled, we can simplify the determinant computation further by introducing energy 
projectors (l7l | and moving to momentum space using Fourier transformation, whereby we find 



lnZ 1>2 = Tr^-^^H^^o-^XJ 2 )/ 3 '^-^^) 2 ))], 

a p ,p 

lnZ 3 = Tr^^H^^-^+J 2 )/? 2 ^ 2 -^) 2 ))] 



(G) 



E 



^E Pia ±(x, E Pta = Jpl + p= + AP , p1 =pi+ P y, else p^ ±a = 2\a\eBn . (7) 



The sum over po = iuj^ denotes the discrete sum over the Matsubara frequencies, n labels the Landau levels in the 
magnetic field which is taken in the z direction and a=0,±l/2 are the possible Q charges of the quarks. 



III. GAP EQUATIONS AND SOLUTION 



Using the following identity we can perform the discrete summation over the Matsubara frequencies 



5>[/3 2 (p 2 - E 2 )] = (3[E + 2Tln(l + e^ E ) = (3f{E) . 



(8) 



Po 
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Then go over to the 3-momentum continuum using the following replacement 

r ^3p 

>^ — > V / a where V is the thermal volume of the system. (9) 

Finally, the zero-field thermodynamic potential can be expressed as 

n 2 A 2 f 00 r] 3 r> 

n B=0 = + -2 -^[f(E+) + f{E;) + 2f(E+) + 2/(25+)] . (10) 

In presence of a quantizing magnetic field, discrete Landau levels suggest the following replacement 



oo 



3 ^2 a n I dp z where a n = 2 - S n0 , (11) 



o (2irf 8tt 2 



71=0 



where a n is the degeneracy factor of the Landau levels (all levels are doubly degenerate except the zeroth level). 
The thermodynamic potential in presence of a magnetic field is given by 

~ ^5>»/ ^[/(^i) + /(^i) + 2 /(<i) + 2 /( £; A,i)] ( 12 ) 

In either case, we can now solve the gap equations obtained by minimizing the (zero-temperature) thermodynamic 
potential Q. 

dfl dfl 
Chiral gap equation : = , Diquark gap equation : = (13) 

Since the above equations involve integrals that diverge in the ultra-violet region, we must regularize in order to 
obtain physically meaningful results. We choose to regulate these functions using a sharp cut-off (step function in 
|p|), which is common in effective theories such as the NJL model [U, [35j], although one may also employ a smooth 
regulator 0, [23|, [34| without changing the results qualitatively for fields that are not too large 2 . It should be noted 
that all free parameters of the system, viz., the bare quark mass, momentum cut-off and quark-antiquark coupling 
constant, are chosen at B = fitting the pion mass, pion decay constant and constituent quark mass in the vacuum. 
Our choice of cutoff implies 

d 3 p \a\eB n ^ f A ' . . r A'~ 



n3 

n=0 



(2n) d Stt 2 



y^a n dp z where n max = Int ; A = ^ A 2 - 2\a\eBn . (14) 

J -A' L2|a|e5J 



From the fact that p 2 > 0, we can calculate the maximum number of completely occupied Landau levels n max . 
For magnetic field B< 10 18 GeV 2 , n max is of the order of 50 and the discrete summation over Landau levels becomes 
almost continuous. In that case, we recover the results of the zero magnetic field case, described in the next section. 
For fixed values of the free parameters, we were able to solve the chiral and diquark gap equations self-consistcntly, 
for B = as well as B large. Before discussing our numerical results, we note the origin of the interdependence of 
the condensates. The chiral gap equation contains only Gs which is determined by vacuum physics, but also depends 
indirectly on Gp/Gs (a free parameter) through A, which is itself dependent on the constituent M = m + a. Our 
numerical results can be understood as a consequence of this coupling and the fact that a large magnetic field stresses 
the qq pair (same Q charge, opposite spins implies anti-aligncd magnetic moments) while strengthening the qq pair 
(opposite Q charge and opposite spins implies aligned magnetic moments). 



2 For example, a smooth cutoff was employed in [34| to demonstrate the De-Haas Van Alfen oscillations in the gap parameter at very 
large B. 
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FIG. 1: Stable, metastable and unstable branches of Q as determined for the case B=0 and i?=0.75 (left panel) and 5=0.05 
GeV 2 and i?=0.75 (right panel) 



IV. NUMERICAL ANALYSIS 



We first fix the NJL model's free parameters; the bare mass m u ^md~5.5 McV, three-momentum cut-off A=0.64 
GeV and G s =5.32 GeV" 2 . These values are chosen to fit three vacuum quantities in the chiral limit: pion mass 
777.^=134.98 MeV, pion decay constant ^=92.30 MeV and the constituent quark mass M(/x = 0)=330 MeV [H SI- 
Although Fierz transforming one gluon exchange implies R=Gd/Gs=0.75 for N c =3, the underlying interaction at 
moderate density is bound to be more complicated, therefore we choose to vary the strength of the diquark coupling 
channel to investigate the competition between the condensates as a function of the chemical potential. In the following 
we first discuss our results for the case with zero external field and then for non-zero external field. 

Case (a) No magnetic field: As long as R is sufficiently small (we will quantify this below) and explicit breaking 
of chiral symmetry is slight, then for a narrow range of fj, we get three solutions to the gap equations. These three 
solutions correspond to stable, metastable and unstable branches of the system. In Fig [1] we plot the values of £1 
corresponding to the three solutions obtained in the small /x-window. The values of the gaps (M, A) for which f2 
is the lowest correspond to the stable solution at any given density. The transition point or the critical chemical 
potential fx c is the point where the first derivative of £1 (and the gaps) behave discontinuously. The superconducting 
gap abruptly becomes non-zero above /j c while the chiral gap abruptly decreases implying partial restoration of chiral 
symmetry (since explicit breaking remains). The behavior of Q and the gaps give us a clear indication of a first-order 
transition. In Fig[2]we plot the gaps versus (/j — /i c )//j, which enables a convenient comparison, on the same graph, of 
the strength of the transition for different values of R=Gd/Gs- The actual value of this critical chemical potential, 
eg. for i?=0.75 (which qualifies as sufficiently small) is /i c =331.2 MeV. Table [J displays /j c and the magnitude of the 
jump in the order parameters at \x c for different values of R. 

With increasing R, the superconducting gap appears at a smaller fj, and rises smoothly from zero, until it becomes 
discontinuous at a critical chemical potential where the chiral gap changes discontinuously (albeit with a smaller 
jump than before). The broken lines in Fig. [21 for R = 1.05, clearly demonstrate this coupled feature of the two 
condensates, with a first-order chiral transition affecting the superconducting gap (see also Huang et al. [l7| for this 
feature). While not obvious from Fig. [2] due to the shift [i — > (fj, — /ic)/V, the chiral phase transition occurs at a 
lower [i than before. This picture does not change qualitatively until we go above a critical value i?=i? c =l.ll; i.e, 
while R < R c , the discontinuities in m and A remain non-zero but decrease smoothly. At R = R c the metastable 
and unstable regions vanish completely. This qualifies it to be a second order phase transition. Huang et al. (l~7| 
have termed the region where the condensates co-exist and vary continuously, as the mixed broken phase, since both 
chiral and (global) color symmetry are broken here. While it should not be confused with a genuine mixed phase, 
since the free energy admits a unique solution to the gap equations in this regime, the width of this overlap region 
increases with increasing R. Above R = R c = 1.11, a and A are smoothly varying resulting in a smooth crossover. 
However there is always a pseudo-transition point (/i cp ) around which fluctuations/variations of both the condensates 
are sharply peaked. The width of these peaks broaden with further increase of R and fi cp moves towards the left with 
increasing R. These results for B = agree qualitatively with the results of Huang et al. [l7| with minor quantitative 
differences in /i c ,Mcp at less than 1.5% level. 
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FIG. 2: Chiral gap (left panel) and Superconducting gap (right panel) as a function of shifted and scaled chemical potential 
(1 — /ic/A 1 ) for different values of R—Gd/Gs- 
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FIG. 3: Chiral gap (left panel) and Superconducting gap (right panel) as a function of chemical potential fi for B=0 and 
B ~ 10 18 G. The value of R=G D /G s =l-2 in both cases. 
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Case (b) Large Magnetic field: Our self-consistent method of solving the gap equations is now applied to 
study the competition between the condensates in presence of a strong magnetic field, at maximum of the order of 
10 18 GeV 2 3 . Since the chiral condensate involves quark spinors of opposite spin and same Q-ch&ige, it is stressed by 
large magnetic fields. The diquark condensate, with opposite spin and Q-charge, is strengthened. Thus, we expect 
and indeed observe, as seen from Fig. [H a strengthening of the competition between the two condensates, resulting 



TABLE I: Critical chemical potential [i c , jumps in the chiral (<S CT ) and superconducting (5a) order parameters at /x c for zero 
and large B-field for various values of R=Gd/Gs- The nature of the transition is also indicated. 



R 


eB=0 


eB=0.05GeV 2 


5 CT (MeV) 


8 A (MeV) 


Mc(MeV) 


Nature 


5a (MeV) 


5 A (MeV) 


Mc(MeV) 


Nature 


0.75 


237.9 


80.6 


331.2 


First order 


158.4 


53.4 


329.5 


First order 


1.05 


115.1 


51.2 


295.3 


First order 


105.6 


66.2 


291.5 


First order 


1.11 








284.9 


Second order 


67.1 


44.4 


280.7 


First order 


1.15 


Smooth 


Smooth 


279.4 


Crossover 


57.0 


32.8 


273.3 


First order 


1.20 


Smooth 


Smooth 


269.7 


Crossover 


40.7 


25.5 


263.2 


First order 



3 We have ignored magnetic catalysis of chiral symmetry breaking, important for extremely large fields B ~ 10 20 G [37| 
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in a qualitative change from the B = case. For moderate values of R < 1.0, the sharp first order jumps in the B=0 
case are weakened and the metastable region shrinks. The effect is more dramatic in the crossover region (large R), 
where the smooth behavior is replaced by a first-order transition (sec Fig [3]). With our choice R — 1.2 (employed 
in Fig. [3]), we find that a smooth crossover in the B = case at fi cp ~ 270 MeV gives way to a discontinuity in the 
chiral and superconducting gaps at \l c ~ 263 MeV for eB = 0.05 GeV 2 (B = 8.5 x 10 18 G). There is only a very small 
change in the transition density for sup erconductivity and (partial) chiral restoration (see Table HI . However, the 
mixed broken phase of Huang et al. [17| is disrupted. The simultaneous appearance of the discontinuity in the chiral 
and superconducting gap for large magnetic field case, at almost the same /z where the condensates have their most 
rapid variation in the B=0 case, is a physical feature and is also cutoff insensitive. Finally, we have checked that 
magnetic fields eB < 0.01 GeV 2 (B < 10 18 G) do not notably alter the competition between the condensates from the 
zero magnetic field case. 



V. CONCLUSIONS 



We have studied the effect of a large magnetic field on the chiral and diquark condensate in a two-flavor super- 
conductor using the NJL model. We have implemented a self-consistent scheme to determine the condensates, by 
numerically iterating the coupled (integral) equations for the chiral and superconducting gap. We have obtained 
results for the nature of the competition between these condensates in two cases: (a) at zero magnetic field, where our 
results are qualitatively the same as those of Ref. [ItJ , and small quantitative differences arise most likely due to the 
choice of free parameters and a different numerical scheme; and (b) at large magnetic fields, where we find that the 
competition increases strongly, causing a discontinuity in the gaps, and disrupting the "mixed broken phase". This is 
a result of the modified free energy of the quarks in the condensate when subjected to a magnetic field. For magnetic 
fields as large as B ~ 10 18 G, the anti-aligncd magnetic moments of the quarks in the chiral condensate change the 
smooth crossover of the chiral transition to a sharp first order transition. The diquark gap also becomes discontinuous 
at this point. For magnetic fields B < 10 18 G, there is no significant effect of the magnetic field on the competition 
between the condensates and zero-field results apply. 

Although we have used an effective model with a sharp cutoff and omitted certain important considerations like 
the finite mass of the strange quark and neutrality constraints on the matter which can induce additional condensates 
that enter the competition, it is interesting to speculate on possible physical consequences of the transition induced 
by the large magnetic field, since it is a strong and probably general feature of competing condensates. At some fixed 
large value of the local field B and in the small density window of the metastable region, cooling of two-flavor quark 
matter below the superconducting T c can result in the formation of domains or nuggets of superconducting regions 
with different values for the gap, even at constant pressure (a detailed analysis requires assessment of surface and 
screening effects in superconducting quark matter [38l. |39|). The chiral gap (or constituent mass) in this mixed phase 
shows the same feature, hearkening back to the DCC idea If the system is formed in a metastable state and 

then drops to the true minimum of the free energy, nucleation of chirally restored and superconducting droplets can 
happen simultaneously. This is a heterogeneous nucleation process since superconducting domains can already exist. 
We can also imagine that external magnetic fields of order eB/^i 2 ~ 1 show some local variation on the microscopic 
scale in the initial stages of the formation of the dense neutron star. In this case, magnetic domains with different 
magnetization can form. Such kinds of nucleation and domain formation will release latent heat that might be very 
large owing to the large value of the magnetic field, serving as a large internal engine for possible energetic events 
on the surface of the neutron star. Such internal mechanisms arc unlikely to occur in a pure neutron star without 



a quark core [28| and could have important applications for high-energy astrophysics 41 1. In addition, it has been 
recently shown that large magnetic fields in strange quark matter introduce anisotropics in bulk viscosities, which can 
change the stability region of the r-mode in such stars [42j . A similar effect may arise in magnetized superconducting 
two-flavor or 2SC+S quark matter, another important link between the microscopic and astrophysical consequences 
of a large magnetic field in quark matter in the core of neutron stars. 
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